Coupling of Low Energy Electrons in Optimally Doped 
E^S^CaCu^Os+tf to an Optical Phonon Mode 

J.D. Rameau, H.-B. Yang, G.D. Gu, and P.D. Johnson 

Condensed Matter Physics and Materials Science Department 
Brookhaven National Laboratory Upton, NY 11973 
(Dated: February 9, 2010) 

Abstract 

Laser based photoemission with photons of energy 6 eV is used to examine the fine details of 
the very low energy electron dispersion and associated dynamics in the nodal region of optimally 
doped Bi2212. A "kink" in the dispersion in the immediate vicinity of the Fermi energy is associated 
with scattering from an optical phonon previously identified in Raman studies. The identification 
of this phonon as the appropriate mode is confirmed by comparing the scattering rates observed 
experimentally with the results of calculated scattering rates based on the properties of the phonon 
mode. 
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I. INTRODUCTION 



In the last decade, with increased energy and momentum resolution, angle resolved pho- 
toemission spectroscopy (ARPES) has emerged as a powerful probe of the mass renormaliza- 
tions, or "kinks", in the electron spectral function associated with the coupling of electrons 
to collective excitations in condensed matter systems.- In particular in many systems it has 
been possible to clearly identify the coupling to phonons, spin waves and charge density 
waves.- In the high T c cuprates on the other hand the identification of the modes responsi- 
ble for the observed mass renormalizations has often proved controversial, particularly when 
different modes occur at the same energy. Experimentally a mass renormalization and asso- 
ciated kink in the dispersion at 70 meV- has been attributed to either spin excitations 3 or to 
phonons 4 . However a recent theoretical study 5 has suggested that coupling to phonons in the 
70 meV energy range and above is, by itself, too small to produce the observed mass renor- 
malization. A mass renormalization observed at 350 meV has also attracted considerable 
attention but again that exact region remains controversial. In the present communication, 
using laser based photoemission~ we present evidence for an additional kink in the nodal 
region of optimally doped Bi 2 Sr2CaCu 2 8+ 5 (Bi2212,T c =91 K). We demonstrate that this 
new kink at low energies exhibits a unique dispersion with respect to the superconducting 
gap in the nodal region of the Brillouin zone. That is, it always appears 8 meV below the 
superconducting gap regardless of the gap's magnitude. Because of this behavior we can 
unequivocally associate an observed mass renormalization in a high T c cuprate with a single 
phonon mode. 

II. EXPERIMENTAL METHODOLOGY 

The ARPES study was carried out with a laser based 6 eV photon source produced as the 
fourth harmonic of a mode-locked Ti: Sapphire oscillator. The oscillator emitted transform 
limited pulses of 2 ps duration at a repetition rate of 105 MHz. The spectral width of the 
fourth harmonic beam was confirmed to be less than or equal to 2.5 meV FWHM with a 
deep ultraviolet imaging spectrometer at the limit of its resolving power. The 6 eV photon 
beam was focused onto the sample with a 100 mm focal length VUV grade CaF 2 lens. The 
sample, a high quality, single crystal of optimally doped Bi2212, was cleaved in situ at a 
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FIG. 1. (Color online) False color intensity maps of the raw (panels a-c) and deconvolved (panels 
d-f) ARPES spectra. The color scaling is depicted at the top of panel c. The intensity scales of 
the raw spectra are normalized to those of their deconvolved counterparts. 

temperature of 10 K. All data were acquired at the same 10 K within five hours of exposing 
the fresh crystal surface. The chemical potential was referenced to a pure gold sample in 
electrical contact with the Bi2212 crystal. Photoemission spectra were recorded with a 
Scienta SES-2002 spectrometer with an overall energy resolution in the present experiment 
of 13 meV and a momentum resolution of .003 A -1 at the photon energies used. 

ARPES spectra acquired at three locations along the Fermi surface are shown in Figure 
[]Ja-c). The location of these cuts in the Brillouin zone is shown in the inset of Figure^. It 
has been shown previously that at a photon energy of 6 eV the antibonding band reflecting 
the bilayer splitting is entirely absent from the spectral. As discussed below the application 
of the Lucy-Richardson method (LRM)^ reduces the broadening due to the instrumental 
resolution, thereby sharpening the spectra considerably. This is illustrated in figures [T] (d-f). 
Following such an analysis the effective energy resolution in the experiment was of the order 
of 5 meV and more importantly for these studies the effective momentum resolution was of 
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FIG. 2. (Color online) Panels a-c show scattering rates r^(w) extracted directly from EDC's of 
the spectra shown in Figure QJl-f, respectively. The dotted lines denote the gap energy relative to 
the midpoint step energy. Red lines are fits to Eq. [T]with a 2 =1.5, 1.1 and 1.2 including resolution 
broadening, impurity scattering and the oJkink = Ar+ ^o dependence of the kink energy. The inset 
of panel a diagrams the locations of cuts in the Brillouin zone. The gap values shown in panels a 
and c have been corrected for the finite energy resolution and kinematic effects, respectively. 

the order of 0.001 A -1 . 



III. EXPERIMENTAL ANALYSIS: THE LUCY-RICHARDSON DECONVOLU- 
TION 

As has been pointed out previously^ ARPES spectra acquired with imaging analyzers 
such as the Scienta 2002 used in the present experiment suffer from a non-trivial finite 
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experimental energy and angular resolution. The finite resolution of an ARPES experiment 
manifests itself in the fact that the spectra acquired, e.g. Figure [T^-c, are an intensity 
map I(k,co) oc R(k,u\k',uj') <8> A(k', u')f{u') where A is the underlying spectral function 
describing the distribution of electronic states in momentum and energy in a sample, f(oo') 
is the Fermi-Dirac distribution for a given temperature and R is a kernel describing the 
resolution broadening in momentum and energy that obscures the spectral function. This 
convolution is responsible for generating the broadening observed in a real ARPES spectrum, 
/. Here the proportionality denotes transition matrix element effects unimportant to the 
analysis. R is taken to be a two dimensional Gaussian distribution with full widths at half 
maximum parametrizing the energy and momentum resolution, respectively. 

While a full description of the how the LRM as applied to ARPES will be presented 
elsewhere^, it suffices here to recall that it is an iterative, essentially statistical "fitting" 
procedure that requires as it's only input the experimentally acquired spectrum I and the 
resolution kernel R which is itself routinely acquired in the course of an experiment. Thus, 
unlike the Maximum Entropy Method (MEM) of spectral deconvolution no a priori assump- 
tions as to the nature of the underlying spectral function are required to extract the spectral 
function. This is highly beneficial when studying phenomena with ARPES for which no 
underlying theory is agreed to exist. When applied to ARPES spectra the LRM produces a 
deconvolved image A that is, statistically speaking, the most likely physical spectral function 
to exist in nature that results in the measured image I when convolved with the inevitable 
instrumental broadening R. 

It is well known that the ill effects of finite resolution on ARPES spectra are most predom- 
inant at and around the Fermi levet^. Indeed, the LRM proves most valuable when called 
upon to correct the measured peak widths of EDO's and MDC's close to the Fermi level. 
This correction has the added bonus of correcting the measured Fermi velocity vp around 
the Fermi level, which tends to be radically increased in an unphysical way by resolution 
broadening. Unfortunately, the widths and dispersions of very low energy electrons, pre- 
cisely where this problem is most acute, are also often of the most fundamental importance 
when comparing ARPES experiments to theory. 

As an example of how the problem of finite resolution can be partially mitigated by the 
LRM in an ARPES experiment we plot in Figure [3] raw and deconvolved EDCs and MDC's 
taken at kp and Ep, respectively, from FiguresfT^ anddH. The sharpening of both the EDCs 
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FIG. 3. (Color online) a) EDC's at kp from the original data of Figure QJ, (open circles) and 
deconvolved data of Figure [Qi (dots), b) MDC's at Ep from the original data of Figure (circles) 
and deconvolved data of Figure [th 1 (dots). 

and MDC's is clearly evident in the comparison. Further, the shift in the peak position of 
the MDC with deconvolution is indicative of the Fermi velocity also being corrected by this 
procedure. The widths obtained from Lorentzian fits of the deconvolved EDC and MDC 
are shown in Figure HI The improvement of both the overall shape as well as widths of the 
peaks is self-evident. 



IV. ANALYSIS OF DECONVOLVED DATA 



As will be shown below, the broadening of the band at binding energies greater than the 
kink energy is in fact the result of an interaction between electrons in the bonding band 
with a well defined low energy bosonic mode. This interaction produces a renormalization 
of the real and imaginary parts of the electron self energy, a " kink" , that itself disperses at 
fixed energy with respect to the bottom of the d-wave superconducting gap. 

Because the observed mode sets in at such low energies and causes rapid changes in 
the band dispersion over a very small energy range an analysis of momentum distribution 
curves (MDC's) does not suffice to extract the imaginary part of the self energy MDC 
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FIG. 4. (Color online) a) Lorentzian fit (black line) to deconvolved data of Figure Eh including 
linear background and Fermi cutoff, b) Lorentzian fit with linear background to deconvolved data 
of Figure 3b. 

analyses were further complicated by the Bogoliubov dispersion acquired by the band with 
the opening of the superconducting gap. Indeed the opening of the gap yields a double 
valued function in momentum up to several meV below the gap edge with a very rapidly 
changing, nonlinear dispersion. The spectra were thus analyzed with EDO's by fitting the 
spectra to Lorentzians with a linear background. While it has been suggested that a more 
complicated, intrinsically asymmetric line shape is required to fit EDC's of the cuprates 
accurately we find that the deconvolved data at low binding energies, exemplified in Figure 
HI are indeed very well fit at both the high and low binding energy side by simple Lorentzians. 
Inverse lifetimes or scattering rates, rV(u;), extracted in this manner are shown in Figure [2] 
as a function of binding energy. At very low binding energies, in the vicinity of the gap edge, 
the lifetime rapidly increases and reaches a maximum at k ~ kp, e.g. where the peak value 
of the superconducing gap is measured. Moving from left to right past kp (in Figure 
[[]) the lifetime decreases while the band dispersion folds back. The increase in peak width 
to the right of kp and folding back of the dispersion reflect the redistribution of spectral 
weight between the occupied and unoccupied branches of the energy symmetric Bogoliubov 
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FIG. 5. (Color online) Gap energies A^ (circles), mode energies (squares) and the difference 
between the two 5^ = $1^ — (diamonds) as a function of Fermi surface angle as defined in the 
inset of Figure [2^. 

quasiparticle dispersion that accompanies the opening of a superconducting gap. A steep 
increase in T around 8.0 meV below the gap edges represents a second pronounced feature 
in the figure. Such a sudden increase in peak width, combined with the kinks visible in the 
spectra in Figure dJi-f, are clear indications that the low energy states are interacting with 
a bosonic mode at a well defined energy. 

We plot the mode energy, Q^, obtained from the theoretical fits described below, the 
superconducing gap, A^, and the energy difference 8% = |fig — A^| as a function of Fermi 
surface angle in Figure [51 Remarkably it is found that over the region of Fermi surface 
investigated in this study 5j: has a nearly constant value of 8 meV. The gap energies, 
denoted by the dashed lines in Figure [2] have in the case of panel (a) been corrected for the 
peak shift caused by residual resolution broadening as indicated from simulation and in the 
case of panel (c) by considering the Fermi surface angle of the kink location as a function 
of Fermi surface angle <f> according to the A(0) = A cos(20) parametrization of the d-wave 
gap with A = 40 meV. Great care must be taken in identifying the correct value of the 
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FIG. 6. (Color online) Schematic representaton of the photohole deexcitation process. At the node 
a photohole (open circle) is annihilated by an electron (filled circle) from the Fermi sea (i) when 
the energy is less than that of the mode energy O. When the photohole resides at a binding energy 
greater than £1 (ii) the decay process involves excitation of a phonon (diagonal line). Away from 
the node, the opening of a superconducing gap A in the single particle density of states shifts the 
absolute value of the binding energies of electrons with sufficient energy for phonon excitation in 
the photohole deexcitation process. The decay process in (iii) shows that electrons with energy 
between A and f] + A do not have sufficient energy to excite the mode. Photoholes at binding 
energies greater than + A (iv) are annihilated by electrons with sufficient energy to excite the 
mode. 

d-wave superconducting gap for spectra acquired off of a high symmetry direction because 
at very low electron kinetic energies the kinematic compression effect^ becomes very severe. 

The fact that the step in the scattering rate is so sharp and that the value of 8^ in 
Figure |5] is approximately constant, even in the presence of an enlarging superconducting 
gap, implies zero momentum transfer or coupling to a q— > bosonic mode^. Physically this 
phenomenon derives from the fact that a photohole injected into the system at a momentum 
k and binding energy u greater than 5r can, in the presence of a zone center mode, only be 
removed by an electron decaying in a vertical, momentum conserving transition. Because in 
the superconducting state all electrons below the gap energy are paired, the lowest energy 
electron that can transition in this manner will always originate from the vicinity of the 
gap edge, in the present case 8 meV away, and so as one scans around the Fermi surface 
the mode energy will always appear at the same binding energy relative to the gap. This 
situation is illustrated schematically in Figure [6j 
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The form and binding energy of St observed in the present ARPES study places a strong 
constraint on the possible bosonic modes involved in the de-excitation. Examination of pub- 
lished IR and Raman spectr a 1516 ' — indicates the presence of a very few, nearly degenerate 
c-axis optical phonon modes in the appropriate energy range, the most likely of which is 
the Raman active Ai g mode previously observed between 58 and 65 cm -1 . In light of the 
theoretical work 1 — and the observed behavior of 5^ we assign this optical phonon mode to 
the 8 meV kink observed in the APRES experiment. This observation is consistent with 
recent temperature dependent work at the node alone^ of optimally doped Bi2212. Clearly 
however a full theoretical understanding of the angle dependent dispersion of this kink re- 
quires more detailed work than was presented in^ as the kink is clearly not pinned to the 
maximum d-wave gap but rather the "local" gap magnitude as a function of Fermi surface 
angle. 

Taking Q = 8 meV and F{oS) a Gaussian of full width at half maximum 1 meV from a 
fit to the Raman data we apply the standard Eliashberg equations^- to calculate the mass 
enhancement factor A and the scattering rate rV(a)) with a 2 as a fitting parameter such that 

A = 2 / y -^du 1 

JO CO 

T(u, T) = 2vr / a 2 F(u') [2n(u/) + f(u' + u>) + f(u' - u)]du' (2) 

<J 

where n(u') and /(u/) are the Bose-Einstein and Fermi-Dirac distributions, respectively. 
The scattering rates calculated in this fashion, shown as fits to the data in Figure [21 agree 
well with the data once shifted to account for the angular dependence of 6?. The fits include 
additional terms reflecting a small, constant impurity scattering component (here set to 5 
meV) as well as the residual energy resolution. No attempt was made in this analysis to 
account self consistently for the presence of the superconducting gap or the presence of the 
kink at 70 meV though we expect these effects to be small, excluding the phenomenological 
shift in binding energy of the kink. From these fits we determine a 2 to be between 1.5 closest 
to the node and 1.15 further away yielding A between 0.52 and 0.4, placing the interaction in 
the weak coupling regime and in good agreement with the average A calculated in reference^. 

As an additional check to the fits of the scattering rates we have applied this kink to 
a simulated low temperature "nodal state". We take the lifetime to be Fermi liquid like 
with T = r e _ e + T Imp + T e _ ph and then calculate the self-consistent real part of the self 
energy by Kramers-Kronig transform. The bare dispersion is set to E(k) = vp(k — hp) with 
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FIG. 7. (Color online) Simulation of 8 meV kink. Relative intensities in arbitrary units are shown 
by the color bar on the right. With the inclusion of even modest resolution broadening mass 
renormalization effects are clearly much more pronounced in the imaginary versus the real part of 
the self energy. 

Vp—2.2 eVA and kp=A6. The impurity scattering in this simulation is set to Tj mp = 5 
meV and the electron-electron lifetime is r e _ e = 2/3((7t/cbT) 2 + uj 2 ) with 2/3=1.5822. r e _ p ^ 
is calculated from Equation [2] with a 2 =1.25 and a Gaussian mode of width 1 meV. The 
quadratic, Fermi liquid form of T e _ e is used for simplicity and at very low energies will not 
deviate substantially from the low temperature cubic form deduced in22. The raw spectrum 
is finally broadened with the energy and momentum parameters from the experiment cited 
above. The result is plotted in Figure [7J Though obviously an oversimplification of the 
true low energy nodal state of Bi2212, the simulation captures the essential features of 
the near Fermi level spectrum we are concerned with. Most importantly, we find that 
the renormalization of the imaginary part of the self energy in the vicinity of the kink is 
much more obvious and dramatic than that of the real part. While the imaginary part of 
the self energy appears as a drastic sharpening of the spectrum upon approaching Ep the 
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corresponding change in the real part of the self energy, even with a limited 5 meV energy 
broadening, is much weaker and will naturally appear as a much more subtle feature of the 
true ARPES spectra in the vicinity of the node. This is precisely what is observed in the 
experiment. It is thus not surprising that previous studies of Bi2212, even at reasonably 
high resolution, have not identified this feature in the spectrum. Certainly the reduced 
momentum resolution of higher energy synchrotron based ARPES studies of Bi2212, in 
addition to the contaminating presence of the bilayer splitting near the node, would make 
this feature otherwise impossible to detect. 

V. DISCUSSION AND CONCLUSIONS 

While the study of electron-phonon coupling in correlated electron materials such as the 
cuprates is of interest in its own right there are several ramifications as regards the more 
general problem of understanding the electronic structure of Bi2212. First, unlike the many 
high energy band renormalizations discussed in the introduction, this finding represents an 
unambiguous identification of such a renormalization with a particular bosonic mode. In 
this case the coupling of electrons in the 2D copper oxide sheets of Bi2212 to an out of 
plane c-axis gerade mode is clearly identified. The identification of this mode from the 
step in the imaginary part of the self energy is supported by the agreement between the 
model calculation of the scattering rate in the presence of such a mode and those measured 
experimentally. 

While the form of 8% shown by this mode is in qualitative agreement with the prediction of 
Ref.— , which takes a much more comprehensive approach to modeling the electron-phonon 
coupling in the presence of a d-wave superconducting gap, e.g. indicating that such a 
coupling is possible a priori, the data bares out an important additional observation. That 
is, a full theory of electron-phonon coupling in Bi2212, in the presence of the d-wave requires 
a proper treatment of the angular dependence of the kink as one scans around the Fermi 
surface. In the reference it was assumed that this effect would be small and one could well 
approximate the shift of the kink energy by pinning it the maximal d-wave superconducting 
gap Ao- Here it is shown that the kink energy close to the node is pinned to the "local" gap 
value A (</>). 

In summary we have reported the observation of coupling to a zero momentum optical 
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phonon. Such coupling can be viewed as the relaxation of the lattice associated with the 
creation of a charged photohole in this ionic material. We note that a similar effect has 
been discussed elsewhere (— and references therein). Such relaxation effects have important 
implications for the transport in such materials. 
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